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近 地 引 力 波 编队 构 型 设计 与 初始 化 方法 研究 


张 继 烨 , 代 洪 华 , 党 朝 辉 ,天 晓 硅 


(西北 工业 大 学 航天 学 院 ,710072 西安 ) 


摘 要 : 近 地 小 型 化 引力 波 探测 编队 研制 周期 短 、 成 本 低 , 适 合 于 引力 波 项 目前 期 的 关键 技术 验证 ， 
具有 重要 的 工程 应 用 价值 。 不 同 于 深 空 环境 , 近 地 引 力 波 编队 受 地 球 非 球形 引力 摄 动 影响 显著 , 星 
间距 离 易 产生 发 散 致 使 任务 失败 。 因 而 相对 于 典型 的 深 空 引 力 波 探测 编队 , 近 地 引 力 波 探测 编队 
在 无 控 状 态 下 的 自然 稳定 编队 构 型 设计 及 其 构 型 精确 初始 化 是 一 个 全 新 的 轨道 设计 与 控制 问题 。 
本 研究 以 二 体 引 力 场 下 的 轨道 动力 学 和 相对 运动 动力 学 为 基础 ,对 正三 角形 编队 构 型 进行 初步 设 
计 。 再 以 构 型 长 期 稳定 为 目标 ,从 摄 动 解析 理论 和 智能 优化 的 数值 搜索 方法 对 摄 动 条 件 下 的 自然 
稳定 构 型 进行 求解 ,为 近 地 引 力 疲 探测 编队 提供 了 稳定 的 构 型 设计 方案 。 最 后 ,针对 探测 构 型 的 精 
确 初始 化 问题 ,探索 了 全 局 最 优 的 微分 校正 策略 ,实现 了 近 地 引 力 波 探 测 编队 的 高 精度 构建 。 
关键 词 : 引力 波 探测 ;正三 角形 编队 ;.J, 不 变性 条 件 ; 数 值 优化 ; 构 型 初始 化 
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了 人 Research on configuration design and initialization method of 
near-earth gravitational wave formation 


ZHANG Jiye,DAI Honghua, DANG Zhaohui ,YUE Xiaokui 
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Abstract: The near-Earth miniaturized gravitational wave detection formation has a short development 
period and low cost. Therefore ,it is suitable for the verification of key technologies in the early stage of 
gravitational wave projects and has important engineering application value. Compared with the deep space 
environment , the inter-satellite distance of the formation of near-Earth gravitational waves tends to fluctuate 
greatly under the influence of the perturbation force, posing a challenge to the natural stable formation 
configuration design and accurate initialization. Based on the orbital dynamics and relative motion 
dynamics under the two-body gravitational field, this paper makes a preliminary design of the equilateral 
triangle formation configuration. Taking the long-term stability of the configuration as the goal,the natural 
stable configuration under perturbation conditions is solved with the perturbation analysis theory and the 
numerical search method of intelligent optimization, which provides a stable configuration design plan for 


the formation of near-Earth gravitational wave detection. Finally, for the precise initialization of the 
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detection configuration, the global optimal differential correction strategy is explored, and the high- 


precision construction of the near-Earth gravitational wave detection formation is realized. 
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optimization ; configuration initialization 


引力 波 是 时 空 曲 率 受 黑洞 合并 等 剧烈 过 程 影 响 
而 产生 的 波动 ,对 人 类 认 知 宇宙 的 形成 和 演化 具有 
重要 意义 。2015 年 9 月 14 日 ,地 面 激 光 干 涉 引 力 波 
天 文 台 (1laser interferometer gravitation wave observato- 
ry) 首次 探测 到 来 自 于 黑洞 合并 所 释放 的 引力 波 信 
号 ,取得 了 引力 波 探 测 的 重大 进展 "。 相 比 于 传统 
地 面 探测 方式 ,天 基 引 力 波 探测 方案 探测 臂 更 长 ,可 
达到 几 十 万 甚至 数 百 万 千 米 ,理论 上 能 够 探测 到 更 
宽频 带 的 引力 波 , 因 而 受到 了 广泛 关注 。 与 地 面 探 
测 原理 相仿 ,3 颗 携 带 迈 克 尔 逊 激光 干涉 仪 的 卫星 
在 空间 中 形成 正三 角形 编队 构 型 ,通过 精确 测量 两 
两 卫星 间 的 光 程 变化 实现 引力 波 探测 ”” 。 近 年 来 ， 
受制 于 任务 周期 与 研发 成 本 ,小 型 化 的 近 地 引 力 波 
探测 项 目 也 被 提出 。 这 类 项 目 研 制 周期 短 成 本 低 、 
适合 于 前 期 的 关键 技术 验证 ,具有 重要 的 应 用 价值 。 
最 具 代表 性 的 如 日 本 的 B-DECIGO 项 目 , 其 编队 尺 
度 仅 为 100 km 并 位 于 近 地 轨 道 高 度 , 主要 用 于 空间 
正三 角形 编队 及 无 拖 忠 控制 等 关键 技术 验证 ,并 能 


算法 对 LISA 自然 轨道 构 型 进行 优化 设计 ,使 编队 星 
间距 离 变 化 在 10 年 内 不 超过 星 间 距离 的 0.99% 。 
文献 [9] 利 用 粒子 群 算法 对 考虑 三 体 引 力 摄 动 、]， 
项 摄 动 的 天 琴 共 轨 星 座 进行 优化 设计 ,使 任务 周期 
内 星 间 臂 长 变化 小 于 1% 。 编 队 构 型 初始 化 方面 ， 
文献 10] 分别 基 于 C-W 方程 设计 了 双 脉 冲 、 四 脉冲 
解析 控制 方法 。 考 虑 真实 模型 的 非 线 性 ,文献 [11] 
利用 微分 校正 修正 末端 位 置 误差 ,使 到 达 精 度 满足 
要 求 。 考 虑 摄 动 影响 ,文献 [12 ] 基于 遗传 算法 和 序 
列 二 次 规划 研究 了 相对 运动 的 轨道 转移 问题 。 文 献 
[13j] 推 导 了 J 摄 劲 和 轨道 面 内 外 耦合 的 轨道 要 素 
偏差 线性 动力 学 方程 ,采用 遗传 算法 和 序列 二 次 规 
划 对 总 速度 增 量 优化 。 但 上 述 方法 均 是 采用 先 优 化 
再 修正 的 思想 ,不 能 保证 获得 全 局 最 优 解 。 

综 上 ,考虑 非 球形 引力 摄 动 下 的 近 地 引 力 波 探 
测 三 星 编队 的 构 型 设计 和 初始 化 问题 ,不 同 于 LISA 
的 以 太阳 为 中 心 引力 场 的 大 天 体 引 力 摄 动 干 扰 问 
题 , 也 不 同 于 天 琴 的 地 、 月 引力 场 共 同 作用 下 的 三 星 


够 辅助 主 项 目 DECIGO 对 引力 波 缘 景 信号 的 精确 探 
测 中 。 此 外 ,考虑 到 探测 精度 ,工程 上 对 正三 角形 编 
队 的 稳定 性 提出 了 要 求 :任务 周期 内 编队 的 臂 长 变 
化 幅度 不 超过 1%'"”。 相 比 于 深 空 环境 , 近 地 空 间 
编队 由 于 轨道 半 长 轴 较 小 ,在 摄 动力 的 影响 下 星 间 
距离 易 呈 现 较 大 的 波动 ,对 无 控 状 态 下 的 自然 稳定 
构 型 设计 提出 了 挑战 。 此 外 ,编队 的 构 型 初始 化 精 
度 也 能 够 影响 编队 稳定 性 ,进而 决定 自然 状态 下 的 


星座 问题 ,是 一 个 全 新 的 轨道 设计 和 控制 问题 。 为 
此 ,本 研究 先 以 二 体 引 力 场 下 的 轨道 动力 学 和 相对 
运动 动力 学 为 基础 ,对 正三 角形 编队 构 型 进行 初步 
设计 。 再 以 构 型 的 长 期 稳定 为 目标 ,从 摄 动 解析 理 
论 和 智能 优化 的 数值 搜索 方法 对 摄 动 条 件 下 的 自然 
稳定 构 型 进行 求解 ,为 近 地 引 力 波 探 测 编队 提供 了 
稳定 的 构 型 设计 方案 。 最 后 ,针对 探测 编队 构 型 初 
始 化 问题 ,本 研究 提出 了 优化 算法 齿 套 微分 校正 的 


星 间 距离 变化 。 为 此 ,本 研究 开展 近 地 引 力 波 空间 
正三 角形 编队 自然 稳定 构 型 设计 方法 与 构 型 精确 初 
始 化 研究 。 

在 编队 构 型 设计 方面 ,文献 [4] 人 研究 了 基于 C-W 
线性 方程 的 编队 构 型 设计 方法 ,给 出 了 空间 圆 直线 
等 构 型 的 详细 设计 流程 。 文 献 [5- 6] 分 别 研究 了 基 
于 常规 轨道 根 数 差 法 和 以 偏 近 点 角 、 平 近 点 角 为 变 
量 下 的 改进 根 数 法 下 的 编队 构 型 设计 方法 。 考 虑 线 
性 模型 精度 ,文献 L7] 基 于 几何 法 设计 了 星 间距 达 
5 x10" km 的 正三 角形 编队 ,以 相对 轨道 面 倾角 为 优 
化 变量 对 构 型 进行 优化 ,设计 精度 可 与 C-W 方程 二 
阶 解 相当 。 考 虑 摄 动 影响 ,文献 [8] 利用 禁忌 搜索 


策略 ,与 前 人 关于 LISA 天 琴 等 任务 的 研究 相 比 ,能 
给 出 多 脉冲 假设 条 件 下 的 全 局 最 优 解 ,并 实现 近 地 
引力 波 探测 编队 的 高 精度 构建 。 


1 正三 角形 编队 初步 设计 


1.1 基于 线性 化 相对 动力 学 的 编队 设计 


引力 波 正三 角形 编队 的 形成 机 理 在 文献 L14] 
中 已 有 详细 论述 。 考 虑 到 编队 尺度 仅 为 百 干 米 并 位 
于 近 地 轨 道 高 度 , 宜 采 用 相对 绕 飞 的 构 型 设计 方式 。 
为 设计 正三 角形 编队 , 先 引 入 一 个 运行 在 圆 轨道 的 
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中 心虚 拟 主星 SC0 ,并 以 SC0 为 圆心 , > 为 半径 作 空 
间 圆 。 空 间 圆 平面 与 SC0 的 轨道 平面 夹 角 为 
60°1"1 ,3 颗 编 队 星 SC1 、SC2 .SC3 以 120° 等 相位 分 
布 在 空间 圆 上 形成 正三 角形 编队 构 型 (图 1)。 以 虚 
拟 主星 为 中 心 建 立 相 对 运动 的 基准 坐标 系 oxyz 。 
编队 卫星 i 相对 于 虚拟 主星 的 位 置 矢量 可 以 表示 为 
pi =Ti -r= [x,y,z] 。 由 于 虚拟 主星 运行 轨道 为 
加 轨道 , 且 主 从 星 距 离 远 小 于 主星 轨道 半 长 轴 。 因 
而 在 坐标 系 oxyz 中 ,从 星相 对 于 虚拟 主星 的 动力 学 
模型 可 近似 用 线性 的 C-W 方程 描述 为 


XT = 2ny + nx +u, 


y = -2nx +u, 和 


区 2 
z=-nz+u, 
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x(t) = xocos(nt) + Fosin( nt) 

y(t) = yocos(nt) — 2xosin( nt) (4) 
2 

2{) = zocos( nt) 十 一 Sin( nt) 
n 


由 从 星 在 运动 过 程 中 始终 处 于 以 虚拟 主星 为 圆 
心 .半径 为 的 空间 圆 上 ,有 


xz (+yY(D +z(D =7 (5) 
由 式 (5) 可 以 获得 补充 初 值 条 件 
Zo = + V3xo, 20 -== (6) 


xo = (r/2)cos0;, Xo =— (rn/2)sing, (7) 
式 中 , 0; 表示 第 i 颗 编 队 星 相对 虚拟 主星 的 相位 角 。 
利用 式 (3) (6) 和 式 (7) 的 初 值 条 件 , 并 使 星相 


其 中 , n = Via 为 主星 绕 地 球 运动 的 平均 角速度 。 
当 不 考虑 控制 力 时 ,方程 的 解析 解 为 


% = 一 (3% + Do]cos(m) + 
n 
0 sin) 十 2 [2x 十 于 | 
n n 


站 po 
y= 2 3% + 2 )sin(nt) Fonly — (2) 
n 


n 

, 2x 
3n 2 十 2) 十 yo 一 一 

n n 


Z 
2 = zocos(nt) 十 sin( nt) 
n 


图 1 正三 角形 编队 示意 图 
Fig.1 Schematic diagram of regular triangle formation 
设 定 相 对 运动 中 心 为 虚拟 主星 , 且 主 从 星 距 离 
不 会 随时 间 发 散 , 有 


1 
-Fm =0 (3) 


基于 式 (3) 的 初 值 条 件 , 式 (2) 可 以 进一步 写 为 


yo + 2nxo = 0, wo 


位 角 相 差 120° , 即 可 实现 正三 角形 编队 的 初步 设计 。 
1.2 能 量 匹配 修正 方法 

1.1 市 空间 正三 角形 构 型 是 基于 一 阶 C-W 方程 
设计 的 ,而 C-W 方 程 存在 的 线性 化 误差 会 导致 主 从 
星 能 量 或 周期 上 的 不 匹配 。 为 了 进一步 缩小 正三 角 
形 构 型 的 稳 态 误差 ,从 能 量 匹配 的 角度 对 上 述 设计 


结果 进行 修正 。 
编队 从 星 的 能 量 可 表示 为 
1 
CAL (8) 
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其 中 ,vj 为 编队 卫星 绝对 速度 矢量 ,其 表达 式 为 


xX-—ny+r, 1 
=ly+Pnxz+r)| = 日 (9 ) 
Z 人 


其 中 :7 为 主星 相对 地 心 和 撩 径 的 大 小 ;7 为 主星 运动 
速度 的 径 向 分 量 。 

由 于 主星 轨道 近 圆 轨 道 , 有 疡 =w,r =0,a 为 主 
星 轨道 半 长 轴 。 将 式 (9) 代 入 式 (8) 可 得 


pA 
= +v2 + ) 六 
= my) +[y tn (rta)’ +2]} (10) 
从 星 能 量 还 可 以 写 为 与 轨道 半 长 轴 相 关 的 形式 
ca = -3 (11) 


将 式 (11) 代 入 式 (10) 中 ,得 
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2 
yo = = KH a = (ny -A n(x+a) (12) 
of Mat +y +7 


式 (12) 即 为 考虑 模型 非 线性 应 用 能 量 匹配 修正 后 y 
向 速度 ,可 进一步 提升 模型 的 精度 。 

分 别 以 五 、L 和 厂 表示 SC1-SC2 .SC1-SC3 和 
SC2-SC3 间 的 距离 。 给 出 虚拟 主星 的 轨道 根 数 如 表 
1 所 示 , 令 3 颗 编队 从 星 初始 相位 9 分 别 为 30°、 
150* 和 270°。 星 间距 离 /= 100 km, 任务 周期 设 定 
为 60 d。 

表 1 虚拟 主星 的 轨道 根 数 


Tab.1 Orbit elements of leader satellite 


虚拟 主星 a/km e MK) 0Q/(°) w/(°) f/(°) 


SCO 42249 1x10-4 20 30 30 0 

根据 式 (3)、(6)、(7)、(12) ,并 利用 相对 运动 
与 绝对 运动 的 转换 关系 '” ,结合 惯性 坐标 系 下 位 
置 .速度 与 椭圆 轨道 根 数 的 转换 模型 “可 获得 每 颗 
编队 卫星 的 轨道 根 数 。 为 使 二 体 情况 下 编队 距离 保 
持 周 期 稳定 ,还 需要 保证 每 颗 卫 星 的 羊 长 轴 相 等 , 即 
aj =w=a。 利 以 上 步骤 获得 的 编队 星 轨 道 根 数 如 
表 2 所 示 。 


表 2 各 颗 编 队 卫 星 的 轨道 根 数 


Tab.2 Orbit elements of follower satellites 


卫星 e A) (2/(°) w/(°) 1(°) 

SC1 6.0x10-4 20.0001 29.802 175.404 214.743 
SC2 7.7x10-4 19.9413 30.100 56.215 333.652 
SC3 6.9x10-4 20.0587 30.099 308.216 81.769 


从 表 2 中 可 以 看 出 ,各 颗 编队 星 的 偏心 率 、 轨 道 
倾角 以 及 升 交点 赤 经 都 较为 接近 。 引 入 纬度 幅 角 的 
概念 ,其 数值 上 等 于 近地点 幅 角 和 真 近 点 角 之 和 , 即 
u = w +/, 解 得 3 颗 编 队 从 星 的 纬度 幅 角 分 别 为 
30. 147° .29. 867°* 和 29.985°, 均 接近 主星 ww =30°。 

进行 60 d 轨道 推演 ,应 用 能 量 匹 配 修正 前 后 的 
星 间 距离 变化 如 图 2、 图 3 所 示 。 

可 以 看 出 编队 星 间 距离 呈现 等 幅 的 周期 性 震 
“ 距离 漂移 幅 值 不 超过 星 间距 离 的 1%。 应 用 修 

条 件 后 , 星 间距 离 波 动 从 原 有 的 0.2% 进一步 缩 
0.12% ,说 明了 能 量 匹配 修正 的 有 效 性 。 
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图 3 采用 外 Ee 下 编队 星 两 两 间距 离 演化 


Fig.3 Formation distance change by applying energy matching 


2 J 摄 动 下 正三 角形 编队 设计 


卫星 编队 位 于 地 球 同步 轨道 高 度 附近 ,这 一 轨 
道 高 度 的 有 项 摄 动 和 月 球 引 力 摄 动 加 速度 大 致 处 
于 同一 量 级 。 但 影响 编队 构 型 的 是 摄 动 引 力 差 , 考 
虑 编队 尺度 以 及 编队 相对 地 球 .月球 的 距离 , 近 地 轨 
道 处 J 摄 动 引 力 差 应 仍 占据 主导 地 位 , 故 和 暂且 只 研 
究 六 摄 动 下 的 近 地 引 力 波 稳定 编队 设计 方案 

卫星 在 J/ 摄 动 下 各 轨道 根 数 的 变化 如 式 
(13) ~(16) 所 示 。 


a =0,6=0,i=0 (13) 


0Q2=- Fn EE (14) 
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SN 
pus 

oO 
避 


-4 人 (6 cosi—-1) (15) 


W = nm Re) (3 oos -1) (16) 


由 于 采用 第 一 章 方 法 所 设计 的 三 星 编队 在 轨道 
倾角 略 有 差异 ,因而 J, 摄 动 影响 下 各 星 升 交 点 赤 
经 、 近 地 点 幅 角 以 及 平 近 点 角 的 变化 速率 并 不 相同 ， 
长 期 将 导致 构 型 发 散 。 


2.1 基于 不 变性 条 件 的 修正 方法 


不 变性 条 件 可 以 抑制 主 从 编队 星 间 距离 发 散 
Sa = zs cos2 -1) le 


i 中 的 3 颗 编 队 卫 星 ， ee J 摄 动 

不 变性 条 件 修 正 求 得 各 颗 卫星 的 轨道 根 数 修 正 量 及 

修正 前 后 距离 偏差 幅 值 的 变化 如 表 3 所 示 , J 摄 动 
不 变性 修正 前 后 星 间 距离 演化 情况 如 图 4 所 示 。 
表 3 利用 不 变性 条 件 的 各 星 修 正 量 与 距离 偏差 


Tab.3 Correction and distance error of follower satellites 


by applying J, invariant condition 


ee 


卫星 8a/m 8i/(°) 距离 偏差 /km 

SC1 0.011 1.88x10- 2.72(relative to SC2) 
SC2 2.725 3.25 x10-4 0.99(relative to SC3 ) 
SC3 -2.750 2.57 x10-4 2.43(relative to SC1) 


105 二 未 修正 世 修 正 后 Wr 


0 10 20 30 40 50 60 


0 10 20 30 40 50 60 
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ud 


图 4 利用 摄 动 不 变 形 条 件 修正 
得 到 的 两 两 星 间 距离 变化 
Fig.4 Formation distance change by applying 


J, invariant condition 
从 上 述 结果 可 以 看 出 ,J 摄 动 不 变性 条 件 对 已 
的 距离 发 散 起 到 了 抑制 作用 但 是 却 扩大 了 的 发 
散 误差 , 且 星 间距 离 漂 移 幅 值 不 能 控制 在 1% 以内。 
原因 是 :(DJ, 摄 动 不 变 性 条 件 是 在 一 阶 近似 条 件 下 
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现象 ,其 基本 原理 是 保证 主 从 星 间 升 交 点 赤 经 相对 
变化 率 为 零 ,并 保证 相同 的 交点 周期 "。 因 而 ,用 
不 变性 条 件 可 以 描述 为 
0 
8T .~ 50+ 8M=0 
若 固 定 主 从 卫星 间 的 相对 偏心 率 ie ,可 以 得 到 
不 变性 条 件 下 相对 轨道 倾角 5i; 和 相对 轨道 半 长 
. 5a 的 近似 表达 式 '" 


(17) 


: 4 
87 = 0 aoti)de (18) 
4e8e 5 3 . Na 
; 引 面 二 十 sin(2D) i (19) 


ey Q@J, 摄 动 不 变 性 需 满 足 式 
nj 会 破坏 原 有 的 空间 圆 构 型 条 件 。 


2.2 数值 优化 方法 


对 于 所 设计 的 正三 角形 编队 ,由 式 (13) ~ (16 ) 
和 升 交 点 赤 经 .近地点 幅 角 和 平 近 点 角 会 发 生长 
漂移 , 且 漂 移 的 大 小 均 只 与 轨道 半 长 轴 、 偏 心率 以 
道 倾 角 有 关 。2. 1 节 的 仿真 结果 也 表明 ,修正 
轨道 半 长 轴 a 和 轨道 倾角 i 能够 改善 星 间 距离 发 散 
现象 。 此 外 , 式 (7) 中 的 相位 角 4% 也 直接 决定 各 星 
的 初始 轨道 根 数 ,因而 也 作为 一 个 优化 变量 。 
选择 各 从 星 a;(io) 、ii(to) 和 bp) 作为 待 优 

化 的 参数 。 考 虑 到 优化 效率 ,本 研究 的 优化 策略 是 : 
先 优化 相位 角 4 ,而 后 直接 对 半 长 轴 ai(io) + 
8ci(z) 和 轨道 倾角 2) + 5i(to) 的 修正 小 量 
ci(p) 和 6i(io) 进行 优化 ,其 中 下 标 i 代 表 3 颗 编 
队 星 序号 。 设 积分 时 间 贡 点 为 shi,… ,bi ,以 L; 


y 


表示 星 间距 离 在 积分 节点 7 处 的 数值 大 小 , M; 表示 


星 间 距离 的 平均 值 , 了 表示 积分 节点 星 间 距离 的 广 
关 清 
ML = 各 (20) 


V = > > es (21) 
基于 以 上 目标 函数 对 SC1 .SC2 和 SC3 的 初始 
瞬时 轨道 半 长 轴 修 正 小 量 、 轨 道 倾角 修正 小 量 和 初 
始 相 位 角 进 行 寻 优 计算 。 相 位 角 9 优化 区 间 设 置 为 
[0°,360°] , 半 长 轴 、 轨 道 倾角 修正 小 量 分 别 设 置 为 
[ -10m,10 mj] 、[ -0.01°,0.01°] ,编队 运行 周期 设 
定 为 60 d。 优 化 目标 为 积分 节点 处 星 间距 离 方差 和 
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最 小 。 利 用 遗传 算法 进行 寻 优 获得 的 各 颗 编 队 星 的 


修正 量 以 及 修正 后 的 最 大 星 间 偏差 幅 值 ,如 表 4 所 
示 。 任 务 周期 内 的 星 间 距离 演化 情况 如 图 5 所 示 。 
表 4 优化 后 的 各 颗 卫 星 轨 道 要 素 修正 量 


Tab.4 Correction and distance error of follower 


satellites by applying numerical optimization 


是 ”00 Sa/m 8i/(°) 距离 偏差 /km 
SCl 160.6 4.25 1.05 x10-3 0.862(relative to SC2) 
SC2 280.6 -2.67 1.16x10-3 0.362(relative to SC3) 
SC3 40.6 -2.07 2.59 x10-3 0.627(relative to SC1 ) 
| 
0 10 20 30 40 50 60 
ud 
100.6 
II ll 
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wd 
时 II 
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wd 
图 5 利用 数值 优化 修正 后 的 两 两 星 间 距离 变化 
Fig.5 Formation distance change by applying 


numerical optimization 

表 4、 图 5 的 结果 说 明了 采用 数值 优化 方法 设计 
正三 角形 编队 的 有 效 性 。 图 5 中 3 颗 编 队 卫 星 两 两 
间距 离 并 不 呈现 相同 的 变化 趋势 ,尽管 在 任务 周期 
内 满足 要 求 ,但 是 长 期 仍然 将 呈现 发 散 特 点 。 其 原 
因 是 :采用 相对 轨道 设计 方式 时 编队 星 的 轨道 倾角 
必然 不 相等 ,由 式 (14) ~ (16 ) , 凡 摄 动 下 的 升 交 点 
赤 经 .近地点 幅 角 以 及 平 近 点 角 的 变化 速率 仍 将 不 
同 ,因而 所 设计 的 正三 角形 编队 长 期 仍 将 呈现 距离 
发 散 的 特点 。 也 侧面 说 明了 ,基于 相对 轨道 动力 学 
方式 所 设计 的 编队 ,在 摄 动情 况 下 往往 需要 后 续 轨 
道 保持 控制 以 维持 相对 稳定 构 型 。 


3 探测 编队 构 型 精确 初始 化 
以 多 脉冲 的 方式 实现 编队 的 构 型 精确 初始 化 。 


在 控制 力作 用 下 ,编队 星相 对 虚拟 主星 的 运动 可 以 
写 为 状态 空间 的 表达 形式 , 即 


=AT+Bu (22) 
由 于 C-W 方程 有 解析 解 ,可 利用 式 (2) ,将 上 述 


方程 进一步 改写 为 状态 转移 矩阵 的 形式 , 即 
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X(1) = Bt ) Kt,) +| B(1,7) Budr 
起 中 : X(1) = [sysz,i, 引 "为 从 是 相对 脱 拟 主星 
的 相对 运动 状态 ; | B(1,7) Budr 表示 控制 力 对 ， 


时 刻 状态 的 作用 。 由 于 脉冲 的 时 间 极 短 , 可 近似 简化 
udt = Ar ,获得 脉冲 作用 力 下 的 编队 星 运动 模型 为 


X(t) = DB(t,to X(to) + bE CERT (24) 


(23) 


( ) 
其 | 中 t,t, 了 区 
D, D,, 


At =t -1t,, B =|0,,3,1;.3] 
其 中 , $ 和 矩阵 的 具体 形式 可 见 文献 [18 ] 。 
当 给 定 末 端 状态 (iy) 和 初始 状态 (io) ， 
式 (24) 可 以 进一步 写 为 
FAV = AX 


| ,mso = B10)B, 


(25) 
式 中 
F = [®t,t) ,GD, (t,t,)], 
AX = [X(i) - Bs,to) X(to)1, 
AV = [An Ay ] 
当 给 定 脉冲 的 施加 时 刻 i 时, 和 矩 阵 为 一 常 值 
阵 ,由 式 (25) 和 广义 逆 理 论 可 知 ,车 n = 2 时 ,满足 
rank(F) = rank(F| AX) =6 , 若 FF- 存在 ,方程 组 
存在 唯一 解 AV = FAKX; 当 ”> 2 时 , 若 满 足 
rank(F) = rank(F| AX) = 6 , 则 方程 组 有 多 解 , 其 
通 解 形式 为 
AV=F*AX+(I-F’F)é 


(26) 


式 中 : PF’ 是 下 的 广义 道 矩阵 渤 为 3 x 1 的 任意 向 
量 。 此 时 , 取 最 小 范 数 解 " 
AV = F':AX (27) 


上 述 求解 过 程 是 在 线性 化 模型 下 进行 的 ,而 实 
际 的 动力 学 模型 是 包含 摄 动 项 的 非 线 性 模型 。 因 而 
需要 对 模型 进行 修正 。 仅 考虑 二 体 非 线性 项 时 ,可 
直接 在 主星 相对 坐标 系 下 进行 微分 校正 。 考 虑 相对 
坐标 系 下 的 J, 摄 动 引 力 差 项 是 和 主 从 星 轨 道 倾 角 、 
纬度 幅 角 相关 的 时 变量 ,积分 求解 比较 困难 ,因而 将 
ee 性 坐标 系 下 进行 。 
惯性 坐标 系 下 ,考虑 外 界 摄 动力 p 作用 下 的 卫 
星 轨 道 动 力学 方程 为 


=- 


7 +p (28) 


当主 从 卫星 在 惯性 系 的 初始 位 置 、 速 度 均 为 已 
知 ,可 以 利用 式 采用 数值 积分 方法 获得 任意 时 刻 的 
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主 \ 从 卫星 的 位 置 .速度 ro.g(t) 、 Vo.g(1) 、 rig(t) 、 
vi.s(t) 。 由 于 编队 星 的 脉冲 施加 是 相对 于 以 虚拟 


性 系 下 。 当 编队 从 星 在 上 时刻 相 对 于 虚拟 主星 的 轨 
道 坐标 系 施加 脉冲 增 量 为 5v1,(71) , 则 其 在 惯性 坐 
标 系 下 的 脉冲 增 量 5v1 .s(t) 为 
Sy1.g(1) = A ,gdv1 1(1) (29) 
其 中 , A, ,s(t) 为 1 时 刻 主星 轨道 坐标 系 到 地 心 惯性 
坐标 系 的 转换 矩阵 ,其 具体 形式 见 文献 [20] 。 
利用 式 (28 ) 和 式 (29) 可 求 得 单 次 速度 脉冲 在 
惯性 系 下 的 表示 。 通 过 对 式 (28 ) 分 段 积分 可 获得 
从 星 末 端 时 刻 的 位 置 和 速度 矢量 六 zt) 、 
vi,s(tj) 。 由 于 构 型 初始 化 的 目标 是 使 从 星相 对 主 
星 达 到 预定 的 相对 位 置 和 速度 zj(tj) 、vi1(ty) ， 


riz(ty) = Ap st) ri et) (30) 
vi(ty) =Agp (1t) [vig(t) —Vo.g(t) 可 
wx|rig(t) —ro.s(ty) | | (31) 


其 中 ,mw 表示 惯性 坐标 系 下 主星 的 角速度 矢量 ,有 
_rig(ty) xvi s(t) 
ri 

通过 以 上 流程 最 终 获得 了 编队 星 末 端 t 时 刻 的 

相对 状态 匀 (tj,AVo) = [riz(sy) ,viz(ty)]，。 则 末 

端 # 时 刻 的 相对 状态 误差 为 

AX(t,) = X(t,) - X(t,, AV,) (33) 

其 中 ,外 (t,) 为 末 段 时 刻 期 望 的 从 星相 对 运动 状态 。 

再 次 利用 式 (27) ,获得 速度 脉冲 的 修正 量 
AVo =F "AX(t,) (34) 
修正 后 的 脉冲 速度 增 量 为 AV, = AV。 + 8AV。 ， 


(32) 
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将 式 (28) ~ (34) 迭 代 循 环 ,直到 获得 符合 末端 误差 
要 求 的 速度 脉冲 AV、 。 

以 上 便 建 立 了 多 脉冲 下 的 精确 编队 初始 化 求解 
模型 。 上 述 过 程 中 ,我 们 假设 阵 给 定 , 但 实际 上 
阵 与 时 间 序 列 {1,} 相关 并 直接 影响 初始 化 的 速度 增 
量 。 因 而 ,将 所 作为 优化 变量 ,优化 目标 为 单 星 转移 
脉冲 速度 增 量 和 最 小 ,优化 表达 式 可 以 写 为 


n 


人 J= > 1Awml 


| (35) 
FAV = AX 

仿真 中 , 设 定编 队 星 和 虚拟 主星 同时 释放 ,编队 
初始 化 时 间 为 86 400 s, 各 颗 编队 星 以 三 脉冲 的 方式 
同时 到 达 预 定位 置 。 各 颗 卫 星相 对 主星 的 末端 位 置 
通过 2. 2 节 数 值 优化 结合 轨道 积分 所 获得 。 利 用 遗 
传 算法 求解 脉冲 施加 的 最 优 时 间 序 列 {1,} ,求解 得 
到 的 脉冲 施加 情况 如 表 5 所 示 , 编 队 从 星相 对 虚拟 
主星 的 转移 轨迹 如 图 6 所 示 。 


s.t. 


SC1 转 移 轨 迹 一 s0 


SC1 目 标点 7 


0 


—60 
Y/km 


图 6 编队 星 转移 轨迹 曲线 


Fig.6 Transfer trajectory of formation satellites 


二 50 Ykm 


表 5 正三 角形 编队 卫星 的 脉冲 施加 情况 


Tab.5 Impulse application of regular triangle satellites formation 


卫星 脉冲 作用 时 刻 /s 单 次 脉冲 增 量 /(m .s …) 总 脉冲 (ms ) 末端 误差 值 

SC1 26 549 [1.036 0.0044 1.707] 4.16 位 置 :1.3 x10 -km 速度 :3.1 x10-?m/s 
68 822 [ -0.5304 0.0046 -0.9289] 
68 861 [ -0.5296 -0.0013 -0.9613] 

SC2 23 [0.5313 -0.2214 0.9707] 4.25 位 置 :2.3 x10- km 速度 :3.2 x10- "m/s 
82 793 [0.7356 0.1160 1.3791] 
82 900 [0.7326 0.1045 1.3683] 

SC3 12 727 [ -1.0938 0.0184 -1.8650] 4.18 位 置 :6.2 x10- km 速度 :6.1x10-2m/s 
55 036 [0.5041 -0.0057 0.8165] 
55 135 [ -0.5016 -0.0145 0.9317] 


从 表 5 中 可 以 看 出 :各 颗 编 队 星 的 脉冲 施加 时 
间 并 不 相同 ,但 末端 位 置 .速度 均 可 达到 较 高 精度 。 
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各 颗 星 转移 到 目标 位 置 所 需 的 速度 总 增 量 仅 4 m/s 
左右 ,可 为 后 续 轨 道 保持 、 重 构 预 留 更 多 燃料 。 


上 由 
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结 论 


研究 了 近 地 引 力 波 编队 构 型 设计 与 精确 初始 化 


问题 ,得 到 以 下 结论 。 


可 以 


1) 编队 设计 方面 ,无 摄 动 下 采用 能 量 匹配 方法 
缩小 稳定 编队 漂移 幅 值 ,使 构 型 漂移 从 原 有 的 


0.2% 降低 到 0. 12% 。J 摄 动 会 导致 编队 构 型 发 
散 , 但 ,不 变性 条 件 不 能 保证 构 型 发 散 得 到 完全 改 


盖 


[= 下 % 


采用 数值 优化 方法 通过 对 各 星 瞬 时 半 长 轴 、 轨 


道 倾角 和 相位 角 寻 优 ,可 设计 出 任务 周期 内 符合 要 
求 的 编队 构 型 ,但 后 续 仍 需要 进行 必要 的 轨道 调整 。 


2) 构 型 初始 化 方面 ,通过 在 智能 数值 优化 算法 中 


嵌 套 


微分 校正 可 以 获得 全 局 脉冲 最 优 解 。 仿 真 结 果 表 


明 , 单 星 仅 需 4m/s 左右 的 速度 增 量 便 可 实现 构 型 高 精 
度 初始 化 ,可 为 后 续 轨 道 保持 、 重 构 预 留 更 多 燃料 。 
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